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We study large-scale, low-energy excitations in the Ising spin glass with Gaussian interactions in 
two-dimensions at zero temperature, using an optimization algorithm to determine exact ground 
states. Periodic boundary conditions are applied. Our results for the fractal dimension of the 
surface, d 3 , and stiffness exponent, 8' , for "droplet" excitations, are in reasonable agreement with 
estimates from "domain wall" calculations, and so support the predictions of the "droplet theory". 
Restricting our analysis to small lattices, we do not find an effective value of 8' close to —0.47 as 
has been recently proposed. The effects of averaging over droplets of different sizes are studied and 
are also found to be too small to give 8' as —0.47 for smaller sizes. Larger corrections to finite-size 
scaling would be needed in three and four dimensions in order for the numerical data used to support 
the "TNT" scenario to be compatible with the droplet theory prediction that the stiffness exponent 
is the same for droplets and domain walls. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q 



I. INTRODUCTION 

Despite extensive studies over many years, there is still 
no consensus on the nature of tb^&pin glass state. Ac- 
cording to the "droplet picture" EtErBo the lowest energy 
excitation, or droplet, of linear size I containing a given 
site has a characteristic energy of order l e where 9 is a 
"stiffness" exponent. Droplets are expected to be com- 
pact but with a surface which has a non-trivial fractal 
dimension d s , less than the space dimension d. For Ising 
spins, the only case considered here, the spin glass state 
persists to finite temperature for d > 3, and in this situ- 
ation one has 9 > 0. On the other hand if the transition 
is at T — 0, as happens in d = 2, then 9 < 0. A com- 
mon way to determine 9 is to compute the characteris- 
tic change in the ground state energy when the bound- 
ary conditions in one direction are changed from peri- 
odic to anti-periodic, which induces a domain wall across 
the system. Although the domain wall, which stretches 
across the system from one side to the other, is differ- 
ent geometrically from a droplet, which closes on itself, 
the droplet picture makes the plausible assumption that 
the energy varies with length scale in the same way for 



calculations 



the two types of excitation. "Domain-wja 
give 9 about 0.20 in tkcee-dimensionsHBlZlH, and about 
0.70 in four dimensipaaa — _In d = 2, a range of domain 
wall calculationsSyy^BBEa give 9 about -0.28, see Ta- 
ble Q, which also shows that the best estimates for d s are 
around 1.27. 

In the alternative scenario, known as "replica sym- 
metry breaking" IRSB) , it is assumed that the 
Parisi theoryBBEJEl of the infinite range Sherrington- 
Kirkpatrick model, applies, at least in part, to short 
range systems. In this picture, the energy of droplets con- 
taining a finite fraction of the system does not increase 
with increasing system size. Although, to our knowl- 
edge, there are no calculations of the domain wall energy 
within the RSB picture, there seems to be no reason why 
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TABLE I: Values of the stiffness exponent 8 and surface 
fractal dimension d s for "domain-wall" calculations in two di- 
mensions. The estimates for 8 are consistently close to —0.28, 
which we shall take to be the correct value (to two decimal 
places) for the purposes of this article. The-^post accurate 
estimates of di-*ire those of Bray and Mooral3, 1.26 ± 0.03, 
and Middletontl 1.27 ± 0.01. 



it should vanish. Hence, it appears that different "9- 
like" exponents are needed for domain walls and droplets 
within the RSB picture. To allow for this possibility we 
will denote by 9 the exponent for the energy of domain 
walls and by 9' the exponent for droplets. Thus, in this 
notation, RSB gives 9' = 0. Furthermore, in RSB, the 
surface of the large-scale, low-energy excitations are ex- 
pected to fill space and so have the fractal dimension 
d s = d. This means that the distributions of both the 
"spin overlap" and "link overlap" , defined in Sec. [H] be- 
low, are "non-trivial", in the sense that they have a fi- 
nite width in the thermodynamic limit, as opposed to the 
droplet picture in which they both become "trivial" , i.e. 
are just delta functions, in this limit. 

Recently, an intermediate scenario, kno wn, . as TNT (for 
"Trivial, Non-Trivial"), has been proposedE3i2j, based on 
numerical work in d = 3 and 4 at zero temperature. 
In this picture one has 9' ~ (< 9), as in RSB, but 
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d s < d, as in the droplet picture. Subsequently, sup- 
port for the TNTj-pjcture was found from finite-T Monte 
Carlo simulationscil. Hacwjever, it has also been criticized 
by Marinari and ParisiEj who argued that d s = d (in. 
which case one obtains RSB). In addition, Middletontll 
has criticized TNT in the opposite sense by arguing that 
averaging over droplets of different sizes, which was car- 
ried out in Refs. ^2] and leads to an incorrect value 
of 9 for the small sizes studied, and that asymptotically 
9' = 9, so the droplet theory is correct. 

Mooreii has also argued that Refs. |22| and ^3] obtained 
the result 9' ~ (< 9), rather than 8 = 9, because cor- 
rections to scaling were neglected. However, the correc- 
tions that Moore considers are different. He argues that 
there is a sub-leading (repulsive) contribution to the en- 
ergy of a droplet, varying as which arises because 
the wall cannot intersect itself. The amplitude of this 
contribution is argued to be larger for a droplet, since 
this has to return to its original location, than for a do- 
main wall which does not have to do this and so would 
have less tendency to intersect itself. Moore argues that 



droplet energy is 
who compute the 



evidence for this contribution to tl 
found in the work of Lamarcq et al.t 
energy of droplets in three dimensions by determining 
the minimum energy excitation above the ground state 
which contains a given number of spins including one 
fixed central spin. Moore further argues that a combina- 
tion of this repulsive piece of the droplet energy which 
decreases with increasing size I, and the leading I term 
which increases with increasing I, can give a droplet en- 
ergy which is roughly constant for the small range of 
droplet sizes studied numerically. By contrast, for do- 
main walls, the amplitude of the repulsive correction is 
argued to be small, and so the correct asymptotic value 
of 9' (i.e. 9) would be observed even for small sizes. 

Further motivation for Moore's ideas comes from re- 
sults in d — 2, where a much larger range of system sizes 
L can be studied than in higher dimensions. While T = 
domain-wall calculations give #._ahout —0.28 as noted 
above, both finite-T simulation£jO and T = studies 
in a magnetic fieldEJ, both of which excite droplets rather 
than domain walls, find 9 about —0.47. A good review 
of this situation is given in Ref. gl| Since 9 < in d = 2, 
both the l e and contributions to the droplet energy 
decrease with increasing L but, Moore argues, could com- 
bine to give an effective exponent of about —0.47 for small 
I crossing over to the asymptotic value of about —0.28 for 
larger sizes. A different perspective on these Jesuits has 
in two dimensions been given by Kawashimao, who ar- 
gues that 9' ~ —0.47 is the asymptotic value for droplets 
and different from 9 obtained £epm domain wall calcu- 
lations. Recently, Picco et al.Ej have studied the first 
excited state, as well as the ground state, in a 2-d spin 
glass, and from their results they infer that 9' ~ —0.46. 

With current algorithms^, it does not seem possible 
to push numerical studies in d = 3 to significantly larger 
sizes than has already been done. However, in d = 2, 
there are more efficient algorithms as we shall see, which 



allow larger sizes to be considered, and the goal of the 
present paper is to investigate large-scale, low-energy ex- 
citations in two-dimensions at T = by the— approach 
of Ref. 23, described below. While MiddletonEZI has also 
applied the method of Ref. ^3] to d = 2, our work differs 
from his because we use periodic boundary conditions, 
as opposed to the "open, with boundary spins fixed un- 
der perturbation" or "link periodic" boundary conditions 
used by Middleton. Since they have no surface, sam- 
ples with periodic boundary conditions are expected to 
have relatively small finite-size corrections. Our choice of 



boundary conditions is discussed more fully in Sec. Ill 



In addition, we also perform an analysis in which we just 
look at the large droplets generated, for each lattice size, 
rather than averaging over droplet sizes. 

Our motivation is to understand in detail finite-size 
scaling for spin glasses in two dimensions. This should 
give better insights into finite-size scaling in higher di- 
mensions, and may help us understand whether the result 
that 9' < 9 found in three and four dimensions is real or 
simply caused by corrections to scaling. More precisely, 
by studying Ising spin glass ground states in d — 2 we 
aim to: 

(i) separate out the effects of averaging over different 
droplet sizes, to see how this influences the values 
of 9 and d s for a system in two-dimensions with 
periodic boundary conditions, and 

(ii) look for the crossover, expected by Mooreii, from 
9 about -0.47 to about -0.28. 

We consider a continuous (Gaussian) distribution of in- 
teractions to avoid the additional complications cmaisg, 
from the degenerate ground state in the ± J modelllJa. 

The model and the technique to generate excitations 
used in Ref. |2^ i s described in Sec. 0. Our choice of 
boundary conditions is discussed in Sec. H . The results 
are presented in Sec. [V and the conclusions summarized 
in Sec. |v} 



II. THE MODEL AND TECHNIQUE 



The Hamiltonian is 



H — JijSiSj 

(id) 



(1) 



where the sum is over nearest neighbor pairs of sites on 
a square lattice, the Si are Ising spins taking values ±1, 
and the Jy are Gaussian variables with zero mean and 
standard deviation unity. The lattice contains N = Lx L 
sites. The techniques presented below work for several 
kinds of boundary conditions, but our results will be 
mainly for periodic boundary conditions. 

We find the large-scale, low energy excitations out of 
the ground state by following the procedure of Ref. |2^. 
First of all we compute the ground state and denote the 
spin configuration by ■ We then add a perturbation 
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to the Hamiltonian designed to increase the energy of the 
ground state relative to the other states, and so possibly 
induce a change in the ground state. This perturbation, 
which depends upon a positive parameter e, changes the 



interactions J, 



by an amount proportional to sf^ Sj°\ 



AH(e) 



'N, 



E* 

(i,3) 



(0) o(0) o c 
3 



(2) 



where Nf, = dN is the number of bonds in the Hamilto- 
nian. In the numerical work we took e = 2. 

The energy of the ground state will thus increase ex- 
actly by an amount AE^ = e. The energy of any 
other state, a say, will increase by the lesser amount 



tween the states 



where q\ 0, °^ is the "link overlap" 
'0" and a, defined by 



(0,a) 



1 



E* 

(id) 



(0) c(0) <?(<*) c( c 
°3 D i °3 



be- 



(3) 



in which the sum is over all the Nj, pairs where there are 
interactions. Note that the total energy of the states is 
changed by an amount of order unity. Fig. [l] shows the 
variation of the energy of the ground state and that of 
an excited state which has an energy 5E above it in the 
absence of the perturbation, as a function of the strength 
of the perturbation. For the value of e shown the excited 
state has become the new ground state. 

We denote the new ground state spin configuration by 
§i°\ and indicate by qi, with no indices, the link-overlap 
between the new and old ground states, i.e. 



<li 



1 o(0) c(0) 5(0) 5(0) 

N °i S j ■ 



Similarly we denote the spin overlap by q, where 



: (o) 5(0) 



(4) 



(5) 



and the change in ground state energy by AE, see Fig. [|. 

From Fig. ([I]), we see that the change in ground state 
energy satisfies the inequalities 



eqi < AE <e. 



(6) 



Furthermore, the unperturbed excitation energy 5E of 
the state which has become the perturbed ground state 
(since 8E + eqi < e, see Fig. pi) satisfies the conditions 



Note that 



< SE < e(l-qi). 
2N vol 



1-9 = 



N 



where N vo \ is the number of flipped spins, and 

2AU-r 

1 ~ qi = ^vT' 



(7) 
(8) 

(9) 



energy 




FIG. 1: A sketch of the variation of the ground state and 
an excited state an energy SE above it, as a function of the 
strength of the perturbation, e. The variation of the ground 
state energy as a function of e is shown by the solid line. 
For the value of the perturbation indicated, the unperturbed 
excited state has become the perturbed ground state. We 
denote by AE the change in ground state energy due to the 
perturbation. 



where -/V sur f is the number of bonds crossed by the domain 
wall bounding the flipped spins, i.e. the surface of the 
domain wall. 

Next, following Ref. we discuss how quantities like 
qi and q are expected to vary with system size. Consider 
first those samples where a large excitation comprising a 
finite fraction of the sites is generated, e.g. those where 
\q\ < 0.5. The average of 1 — q over these samples given 
by 



l. 



(10) 



where [• • -] av denotes an average over samples and the 
prime denotes a restricted average over those samples 
with \q\ < 0.5. From Eq. (J^j we see that the average of 
1 — qi over these samples given by 



n „ V T-{d-d s \ 



(11) 



Such a large excitation costs an energy of order L e . 
This must be more than compensated for by the en- 
ergy gained from the perturbation (relative to the ground 

state), e(l — qi), which, from Eq. ([|), is proportional to 
eL -(d-d 3 )^ Hence 

assuming a distribution of droplet en- 
ergies with a finite weight at the origin, the probability 
that such a large excitation can be created is given by 



P(\q\ <0.5)~L 



-{6' + (d-d s )} 



(12) 



We now consider results of averaging over all samples 
rather than just those with large droplets. Eq. ( ^2[ ) shows 
that there is probability ~ L~{® +( d - d s)} to have 1 — q 
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of order unity, and so, assuming that samples with large 
droplets dominate the average, we have 

[1 - g]av ~ L-^ d - d ^ . (13) 
Similarly, we find 

[l-®U~£-«, (14) 

where 

M =6' + 2{d-d s ). (15) 

Note that that the ratio of the quantities in Eqs. jl^ ) 
and ([l4] ) is independent of 9', i.e. 

R= [ } ~ qi } av ~ L^ d - d '\ (16) 

Eq. (0) shows that the droplet excitation energy SE, 
see Fig~|j], scales in the same way as 1 — qi, i.e. like 
Eq. (|ll]) if only samples with \q\ < 0.5 are studied, and 
like Eq. (|l4|) if all samples are averaged over. 

In general we can extract d — d s from the surface to 
volume ratio of the excited droplets of spins, but to get 
0' we also need information about the probability that a 
large droplet is excited. 

III. CHOICE OF BOUNDARY CONDITIONS 

In those two-dimensional models where the bonds form 
a "planar graph" , which unfortunately excludes peri- 
odic boundary conditions in both directions, very ef- 
ficiest, jWjJypamkJ time "matching" algorithms can be 
uscdllil'OaEaOEa to determine the ground state. As a 
result very large sizes can be studiec£3. By contrast, in 
three or higher dimensions the calculation is knownEZI to 
be NP-hard which means that a polynomial-time algo- 
rithm is unlikely to exist. 

In Fig. H we show the excitation generated, using the 
matching algorithm, for a 200 x 200 sample with peri- 
odic boundary conditions on the vertical edges, and open 
boundary conditions on the horizontal edges. Unfortu- 
nately, the boundary of the resulting excitation is not a 
droplet but a single domain wall which spans the whole 
system. This structure would not have been possible if 
periodic boundary conditions had also been applied along 
the top and bottom edges. We find that a single do- 
main wall is generated for a significant fraction of the 
samples with open-periodic boundary conditions, and so 
these boundary conditions are not appropriate for a study 
of the differences between droplets and domain wallscj. 

One possibility to prevent the occurrence of domain 
walls is to first generate a ground state with open bound- 
ary conditions, p-and then fix the orientations of the 
boundary spinsEJ relative to each other by introducing 
very strong bonds in favor of their relative ground-state 
orientations. This still allows us to use the fast matching 




FIG. 2: An excitation generated for an L — 200 system 
obtained using the matching algorithm. The vertical edges 
have periodic boundary conditions while the horizontal edges 
have open boundary conditions. With the matching algorithm 
it is not possible to have periodic boundary conditions in both 
directions. The line shows the boundary separating the region 
where the spins are flipped when the perturbation is applied 
(denoted by "J.") from the region where they are not flipped 
(denoted by "f"). One sees that, in this case, the excitation 
is a single domain wall rather than a droplet. 

algorithm and hence to treat large system sizes. We have 
only tested this approach briefly, since we found that usu- 
ally only very small droplets are generated, with almost 
none spanning more than one quarter of the system (i.e. 
with \q\ < 0.5). The reason is that by fixing the bound- 
ary spins an effective repulsion of the domain wall from 
the boundary is introduced. This reduces the sizes of the 
droplets. Furthermore, it is also likely that this repulsion 
affects the morphology of the domain walls. 

We therefore decided to apply periodic boundary con- 
ditions in both directions which, unfortunately, means 
that the matching algorithm cannot be usedcil. Instead 
we use the spin_glass ground state server at the Univer- 
sity of Colognecj. This uses a "branch and cut" algo- 
rithm for which, as we shall verify, the CPU time in- 
creases faster than a power of the system size, so we 
cannot study such lpge systems as are possible with the 
matching algorithmtia Nonetheless, the implementation 
of the branch and cut algorithm on the Cologne server is 
very efficient, so we can still study quite a large range of 
sizes, in practice L < 64. 

Fig. H shows a large droplet excitation created for an 
L = 64 sample with periodic boundary conditions. In 
some samples where large-scale excitations are produced, 
the region of flipped spins wraps round the sample ere- 
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FIG. 3: A droplet excitation generated for an L = 64 system 
obtained using the branch and cut algorithm of the Cologne 
spin glass server. Periodic boundary conditions are applied in 
both directions. The dots indicate the lattice sites. Clearly 
the region where the spins are flipped forms a droplet. 
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44 
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2099 


2015 


1763 



TABLE II: The number of samples N S!Lln p for each lattice size 
L. 



ating a "sponge-like" excitatppljo, which is, in effect, 
made up of two domain wallsEl An example is shown in 
Fig. |J. Later we will discuss the role played by sponges, 
finding that their properties are very similar to those of 
large-scale, non-sponge excitations. 



IV. RESULTS 

All the numerical results are for a strength of the per- 
turbation e = 2. Asymptotically, the results should be 
independent of the value of e. However, if we make e 
too large, then the asymptotic behavior will only be seen 
for large sizes, and if e is too small the statistics will 
be poor. Based on experience with other modelsEj, the 
choice e = 2 seems to be a reasonable compromise. In 
Table || we show the number of samples studied for each 
lattice size. 




FIG. 4: A "sponge-like" excitation generated for an L — 64 
system with periodic boundary conditions in both directions. 
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FIG. 5: A log-log plot of the fraction of samples for which the 
spin overlap between the perturbed and unperturbed ground 
states satisfies \q\ < 0.5. According to Eq. dl2| ) the slope 
should be —[9' + (d — d 3 )}. The inset is a log-log plot of the 
fraction of samples for which not only is \q\ < 0.5 but also the 
excitation is not sponge- like, see Fig. (Q). The slope is almost 
the same as that for in the main part of the figure, indicating 
that a roughly fixed fraction, about 55%, of excitations with 
\q\ < 0.5 are not sponges. 
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FIG. 6: A log-log plot of the average of 1 — qi over those 
samples for which the spin overlap between the perturbed 
and unperturbed ground states satisfy \q\ < 0.5. According 
to Eq. ([ll]) the slope should be —(d — d s ). The inset shows 
results averaged over samples where not only is \q\ < 0.5 but 
also the excitation is not sponge-like, see Fig. H. The slope 
is virtually the same as that in the main part of the figure, 
indicating that sponges have similar properties to large non- 
sponge, droplet excitations. 



A. Average over samples with large droplets 

Firstly, we consider results obtained by averaging only 
over those samples where large droplets were created, i.e. 
the spin overlap satisfies |g| < 0.5. Fig. || shows results 
for the probability that this occurs, while Fig. ^| shows 
the average of 1 — qi for those samples. Fitting the data 
according to Eqs. (O) and (III]) we find 



d — d s 
9' + (d-d s ) 



0.715 ±0.003 
0.396 ±0.015, 



which gives 



d s = 1.285 ±0.003 
9' = -0.319 ±0.015. 



(17) 



(18) 



The value of d s agrees very well with best estimates of 
around 1.27, see Table while the value of 9' is a lit- 
tle more negative than the best estimates for 9 from do- 
main wall calculations, which is about —0.28, see Table §. 
However, the difference is quite small, only about two 
standard deviations. Hence we see that our results are 
consistent with the droplet theory being correct in two 
dimensions. 

The insets to Figs. || and ^ show data for those samples 
where the excitation is not only large, i.e. \q\ < 0.5, but 
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Qi 


Q 2 


4 < L < 64 
4 < L < 16 
16 < L < 64 


1.285 ±0.003 
1.299 ±0.005 
1.285 ±0.008 


-0.319 ±0.015 
-0.288 ±0.029 
-0.298 ± 0.043 


0.006 
0.037 
0.405 


0.253 
0.050 
0.639 



TABLE III: Fits for the data for those samples where |g| < 
0.5 for different ranges of sizes, Qv and Q2 are the quality 
factors, as conventionally defmedc3, of the fits for [1 — qi]'^ 
and P(\q\ < 0.5) respectively. Note that Qi is quite small for 
the fit with all sizes. This is because the error bars are small 
and there are clearly some systematic effects. For example, if 
we omit the L=4 point then Qi becomes 0.129, which is much 
better. However, this then gives 9 — —0.347 ± 0.019, which 
is further from the accepted value of about —0.28 but agrees 
within the error bars with the value in the table. 



is also not sponge-like, i.e. is not of the form shown in 
Fig. [|. The fraction of samples with \q\ < 0.5 for which 
the excitation is not sponge-like is roughly independent 
of size at about 55%. This is why the slope in the inset 
to Fig. is very similar to that in the main part of the 
figure. Furthermore, the slope of the data for the average 
of 1 — qi in the inset to Fig. |^ is very similar to the average 
over all samples with \q\ < 0.5 (shown in the main part 
of the figure), indicating that the properties of sponge- 
like excitation are similar to those of large non-sponge 
droplets. 

We have also analyzed the data including just the 
smaller sizes, 4 < L < 16, and just the larger sizes, 
16 



< L < 64. and the results are shown in Table III. 
The values of 9' are, given the error bars, scarcely more 
negative than the best accepted value of —0.28. One 
does not see a systematic trend over this range of sizes. 
In particular we do not see the crossover, expected by 
MooreEj, from a value of about —0.47 at smaller sizes to 
about —0.28 at larger sizes. 

From Eq. ([?]) we expect that the unperturbed excita- 
tion energy SE, of the droplet created by the perturba- 
tion, see Fig. [I], should scale in the same way as 1 — qi, i.e. 
with an exponent d — d s if restrict the average to those 
samples with \q\ < 0.5. The data, plotted in Fig. |^, gives 
d-d s = 0.683 ± 0.011 as compared with 0.715 ± 0.003 
obtained from the data for 1 — qi in Fig. |(| The difference, 
which is quite small, is presumably due to corrections to 
scaling. 

We should make clear why the average of SE induced 
by the perturbation varies as L~ lyd ~ ds ^ (if we only in- 
clude samples with \q\ < 0.5), whereas the characteristic 
droplet energy varies as L e . The point is that there is a 
distribution of droplet energies which is assumed to have 
a finite weight at the origin. Thus, some samples will 
have a large-scale droplet with energy much less than 
L e . In fact, the perturbed ground state will only have a 
system-size droplet for those samples where the droplet 
energy is comparable to the energy which can be gained 
from the perturbation, i.e. L~^ d ~ ds \ The probability 
of this occurring is therefore L~{® ag shown in 
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FIG. 7: A log-log plot of the average of the unperturbed 
excitation energy, 6E, see Fig. |l[ over those samples for which 
the overlap between the perturbed and unperturbed ground 
states satisfies \q\ < 0.5. According to Eqs. (R) and (|lj) the 
slope should be — (d — d s ). 



Fig. |. 

Our expectation in looking at samples with \q\ < 0.5 
is that there is a single large droplet which dominates, 
rather than having several smaller droplets of compara- 
ble size, and we have verified that this indeed the case. 
We have separated out the different droplets for each 
sample and looked at averages over those droplets con- 
taining 7V/4 spins or greater (from Eq. (||) this corre- 
sponds to |<7| < 0.5). As an example, for L = 64, out 
of 1763 samples, 245 of them had \q\ < 0.5 and 241 
had a single droplet which would give \q\ < 0.5 by it- 
self. If we repeat the finite-size scaling analysis for these 
large droplets we find d — d s = —0.715 ± 0.003 and 
9' + (d — d s ) = —0.392 ± 0.015, in very close agreement 
with the results in Eq. (|l7|). 



B. Average over all samples 

Next we discuss the results obtained by averaging over 
all samples, rather than just those with |q| < 0.5. Fig. ^ 
shows results for the ratio of the average of 1 — qi to the 
average of 1 — q, and Fig. || shows results for the average 
of 1 — q. Fitting the data according to Eqs. ( |l6| ) and JT3| ) 
we find 



d — d s 
9' + (d- d s ) 



0.674 ±0.004 
0.320 ±0.010, 



FIG. 8: A log-log plot of the ratio [1 — g;] a v/[l — ?]av against 
L obtained by averaging over all samples. From Eq. (hq) the 
slope should be —(d — d s ). 



which gives 



d 3 = 1.326 ±0.004 
0' = -0.354 ±0.010. 



(20) 



Analyzing the data for different ranges of sizes gives the 
results shown in Table IV . We note that the values for eL 



are larger than those obtained from analyzing just large 
droplets, see Table III, and the values for 9' somewhat 



more negative. This is presumably due to the effects of 
averaging over droplets of different sizes, as discussed by 
MiddletonO. Indeed, our results for jjbi = 2{d — d s ) + 0' , 
see Table |v[ a gree quite well with those for smaller sizes 
in Fig. 4 of Ref. [17], even though our boundary conditions 
are different from his. Nonetheless, the variation of expo- 
nents with range of L is4oo small to yield 6' about —0.47 
as suggested by MooreEa, thought the trend (9' more neg- 
ative for smaller lattice sizes and when we average over 
different droplet sizes) is in the right direction. 

From Table Cvl we find 9' = -0.364 ± 0.021 when 



averaging over all samples and just including the smaller 
sizes. This is about —0.08 ± 0.02 different from the best 
value of —0.28, see Table |. Furthermore, a reanalysisEZI 
of the data of Ref. |23|, just keeping those samples with 
small \q\, i.e. similar to our analysis which yielded the 
fits in Table III , found the samej-discrepancy in three- 
dimensions as the original analysiso which averaged over 
droplets of all sizes. We see that the discrepancy in the 
relative to —0.28, is much smaller 
even if the fits are restricted to 



value for 9' in Table fl] 
than that in Table 



IV 



(19) 



smaller sizes. Hence a much larger effect would be needed 
in three dimensions to explain why the numerics found 
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FIG. 9: A log-log plot of [1 — g] av a gain st L obtained by 
averaging over all samples. From Eq. (113) the slope should 
be -[6' +(d-d 3 )]. 




FIG. 10: A plot of the size AT vo i and surface area Asurf of 
individual droplets for L = 64. A fit gives d s /d — 0.632. 





d s 


8' 




4 < L < 64 
4 < L < 16 
16 < L < 64 


1.326 ± 0.004 
1.340 ± 0.007 
1.322 ±0.011 


-0.354 ±0.010 
-0.364 ±0.021 
-0.317 ±0.030 


0.995 ± 0.008 
0.956 ±0.017 
1.039 ±0.024 



TABLE IV: Fits for the data obtained by averaging over all 
samples for different ranges of sizes. Here /i; = 2(d — d s ) ± 8' . 



9' ~ if the correct result is 9' = 9 ~ 0.20. Io four- 
dimensions there is an even greater discrepancy^, 9' ~ 
0,6» ~ 0.70. 

To get an independent estimate of d s we have also 
looked at the size, N vo \, and surface area, N sux f, of in- 
dividual droplets for a given size. Fig. |l0| shows data for 
the largest size L = 64, and similar results are found for 
other sizes. Doing a least squares fit gives d s /d — 0.632, 
or d s — 1.264, in very good agreement with the best es- 
timates of around 1.27, see Table |[ 



C. Computational Complexity 

Finally we discuss how the CPU time depends on the 
system size. Fig. [y] shows a log-log plot of the aver- 
age CPU-time on the Sun Ultra 10/440 of the Cologne 
ground state server against L. The clearly visible up- 
wards curvature indicates that the time increases faster 
than a power of L, i.e. the algorithm is non-polynomial, 
as expected. The inset to the figure shows that, empir- 



ically, the CPU time varies like exp(const. for the 
larger sizes studied. However, we do not know the true 
asymptotic dependence of the CPU time on system size. 



V. CONCLUSIONS 

We have studied large-scale, low-energy excitations at 
T = in the Ising spin glass in two-dimensions with pe- 
riodic boundary conditions. We have particularly looked 
at the dependence of d s and 9' on whether or not we av- 
erage over droplets of different sizes, and on the range 
of system sizes that are included in the fits. Averaging 
just over those s ampl es in which lar ge droplets are gen- 



erated, see Sec. [IV A| and Table \ni[ the results do not 
vary much if we consider different ranges of system sizes, 
and are consistent with the droplet theory prediction that 
6' = 8. If we average over all samples, see Sec. IV B and 



Table IV, which, for each lattice size, involves averaging 
over a range of droplet sizes, the results vary more as a 
function of the range of system size considered. However, 
the effect is too small to give an effective raJue of 8' close 
to —0.47 as has been recently proposed^, though our 
result in line 2 of Table IV does lie between —0.47 and 
the expected value of about —0.28. It is possible that 
different ways of generating droplets will yield different 
effective values of 9' for small sizes, and this may be the 
reason why Refs. HJ1II] nnd d ' - "°- 47 - 

We have also considered the role of sponge-like exci- 
tations, see Fig. [|, finding that their properties are very 
similar to those of large, non-sponge droplets, see Fig. |j. 

Overall, our results support the droplet theory for the 
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FIG. 11: A log-log plot of the average CPU time per sample 
as a function of the system size L. The upwards curvature 
shows that the CPU time increases faster than a power of L 
as expected. The inset shows that the CPU varies roughly as 
exp(const. VL) for the larger sizes. The dashed line is a guide 
to the eye. 

two-dimensional spin glass and do not provide a^stralght- 
forward way to reconcile the numerical results£3'E2l used 



to support the TNT scenario with the prediction of the 
droplet theory that the stiffness exponent is asymptoti- 
cally the same for droplets and domain walls. For this 
prediction to be true, larger corrections to finite-size scal- 
ing would be needed in three and four dimensions than 
in two dimensions. 
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